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Abstract. This article proposes a way to improve the presentation of histograms where data are compared 
to expectation. Sometimes, it is difficult to judge by eye whether the difference between the bin content 
and the theoretical expectation (provided by either a fitting function or another histogram) is just due 
to statistical fluctuations. More importantly, there could be statistically significant deviations which are 
completely invisible in the plot. We propose to add a small inset at the bottom of the plot, in which the 
statistical significance of the deviation observed in each bin is shown. Even though the numerical routines 
which we developed have only illustration purposes, it comes out that they are based on formulae which 
could be used to perform statistical inference in a proper way. An implementation of our computation is 
available at: http://svn.cern.ch/guest/psde/. 

PACS. 02.70.Rr General statistical methods 



1 Introduction 

Most analyses compare the observed data to the expectation resulting from a theoretical model, like the Standard 
Model (SM), or some other hypothesis, like a best fitting function, or a Monte Carlo simulated distribution. It is 
common to make histograms in logarithmic scale, because their contents span orders of magnitude. Since differences 
are difficult to see in a logarithmic scale, an inset plot is often made at the bottom of the histogram. In this article, 
we focus on this inset and propose a way to make it intuitive and accurate. By "intuitive" we mean that it should 
make it obvious which bins contain an excess of data and which contain a deficit, while significant deviations should 
look more striking than insignificant ones. By "accurate" we mean that it should represent the actual significance of 
the deviation in each bin, rather than some approximation. 

Below we define statistical significance in two probability models, Poisson and binomial. In Section 2 we show a 
sequence of presentation options that could be used (and most of them have been used), explaining their strengths 
and limitations. This will motivate, through a series of incremental improvements, our final proposal, which is given 
in Section 2.4. Finally, if the expectation is only known within some uncertainty, this uncertainty can be taken into 
account in the comparison to the data, as explained in Section 3. In this paper, we use the ROOT framework [1] to 
produce the plots and provide suggestions for implementing our formulae with this widely used software for data 
analysis and visualization, although our recommendations can be easily implemented within other frameworks too. 

1.1 Definition of statistical significance 

The key notion, which is needed to define the statistical significance, is the probability of finding a deviation at least 
as big as the one observed in the data, under the assumption that the chosen theoretical model describes our system. 
This probability is commonly called the "p- value" and it usually spans several orders of magnitude. For this reason, it 
is convenient to translate the p- value into a "z- value" , which is the deviation at the right of the mean of a Gaussian 
distribution, expressed in units of standard deviations, which would correspond to the same p-value. The equation 
which maps each p-value to a z-value, and vice versa, is: 




(1) 
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Fig. 1. Relationship between p-value and z-value. 



wfiicfi is shown in Fig. 1. In ROOT, the z-value can be computed in one line of code, using the inverse error function: 

sqrt(2.) * TMath :: Erf Inverse(l. - 2. * pvalue) 

A z-value > corresponds to a p-value < 0.5, and negative z-values correspond to p-value > 0.5. Significant 
deviations are characterized by quite small p-values, corresponding to z-values > 3. For example, it is common to refer 
to a p- value = 2.87 x 10~^ as a "5cr effect", meaning that the corresponding z-value is 5. Such deviations are usually 
considered very significant and a z-value of 5 or more is conventionally required in High Energy Physics to claim a 
discovery [2]. On the other hand, a z-value which is less than 1-2 units represents a common statistical fluctuation, 
which is uninteresting. The fluctuations are even more probable (i.e. less interesting) when the z-value is negative. 

The p-value depends on the statistical distribution that the data are supposed to follow, the so-called "probability 
model" . The cases of Poisson and binomial distributed data are most common in experimental physics, and are 
addressed in this note. 



1.2 Poisson model 



This is the most common case, where event counts are plotted in each bin and the number of entries in each bin follows 
a Poisson distribution (e.g. Ref. [3]). If S G M events are expected in a bin, the probability of observing D E N events 
is 

D\ 



P{D\B) = Poi(i:'|B) = — e 



[B is also the variance of D). The Poisson p- value is 



E 



D-l 



-B 



-B 



p- value — < 



_ ) n=D 



E 

n=0 



B"" 



7 I 



D> B 



D < B 



(2) 



The above sums are simplified thanks to the identity 
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where 



r{D,B) = 



i^-^e-* dt. 



(4) 
(5) 
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The ratio between the upper incomplete Gamma function r{s,x) and the Gamma function r{s) in the previous 
equation is known as the upper regularized Gamma function 

Q{s,x) = ^^j^ = l-P{s,x) (6) 

r{s) 

where P(s, x) is the cumulative distribution function for Gamma random variables with shape parameter s and scale 
parameter 1 [4]. In ROOT this function is available as 

Q{s,x) = ROOT: :Math: : inc_gamina_c (s ,x) , 

such that 

[ 1 - (5(1?, S) = ROOT: : Math: :inc_gainma_c(D,B) ,D>B 
~ 1 + 1,5) = ROOT: :Math: :inc_gainma_c(D+l,B) ,D<B 

One should notice that it is not always true that the bin population follows the Poisson distribution. This does 
not happen when the total number of entries in a histogram is not a random variable, in which case the events in 
each bin are distributed accordingly to the multinomial distribution, which means that the bins cannot be considered 
statistically independent as we assumed in this section. In the extreme case of a fixed total number of entries and 
only two bins, the content of each bin follows the binomial distribution: if there is an excess of counts in one bin, with 
respect to the expectation, then the other bin must have a deficit in order to preserve the total number of events. 



1.3 Binomial model 

If the plotted quantity is a Bernoulli success rate, then it follows a binomial distribution. A common example would 
be a trigger efficiency "turn-on curve" like Fig. 1 in Ref. [6] , where it would be informative to show how significantly 
the observed trigger rate differs from the simulation. Another example is F^^nijj) in Ref. [5]: if this event rate were 
significantly different than the SM prediction, this could be a sign of certain new physics. 
Let us define the following notation: 

n (z N: the number of events initially observed in a bin, a subset of which will survive the selection; 
fcobs € {0, 1, ■ • . , n}: the observed number of selected events; 

e S K: the expected success rate, i.e. the selection efficiency, to which we compare the observed success rate kohs/n. 
The binomial distribution gives the probability of observing k events passing the selection out of the initial n events: 



P(fc|n, e) = Bi(A:|n, e) = i^^j e'^(l - e)""* (8) 

The expected number of surviving events is E[k\n, e] — ne, with variance y[fc|n,£] — ne{l — e). Finally, the Binomial 
p-value is defined as 



p- value = < 



> ne 



' n 



ri=0 



The cumulative distribution function of the binomial model can be represented in terms of the regularized incom- 
plete Beta function [4] Ix{a,b) — Jq ^^^^{1 — t)''^^ dt/B{a,b) where the Euler's Beta function can be represented in 
terms of the Gamma function as B{a,b) = r{a)r{b)/ r{a + b). The regularized incomplete beta function is available 
in ROOT as 

Ixia,b) = ROOT: :Math: : inc_beta(x, a,b) , 

which makes it possible to implement what follows. 

Let us start with the deficit case fcobs < ne, for which the p-value coincides with the cumulative distribution 
function: 

p-value = P{k < fcobs) = h-e{n - fcobs, fcobs + 1) , fcobs < ne . (10) 

For the case of an excess, we use two identities. The first is Ix{a, &) = 1 ^ Ii-x{b, a) and the second comes directly from 
the binomial distribution and is the starting point of the following chain of equations, whereas the previous identity 
is used in the last step: 



P(fc > fcobs + 1) = 1 - P(fc < fcobs) = 1 - Il-e{n - fcobs, fcobs + 1) = 4(fcobs + l,n~ fcobs) 
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Table 1. Observed events (D) and expected events (B), used for demonstration. 



from which we immediately find the p-value for the excess: 

p-valuc = P(fc > kohs) = 4(fcobs, - ^obs - 1) , ^obs > ne (11) 
which is defined for fcobs < n (in case kohs — n one has P{k > n) — P{k — n) — Bi(n|n, e) — e"). 

1.4 Unfolded data 

If the histogram is the product of some unfolding procedure, things are much less clear, because the bin populations do 
not follow a well-known and generally applicable probability distribution. In addition, the correlations across adjacent 
bins are unavoidable. It is unclear even if one should try to compute, bin by bin, the significance of the difference 
between unfolded data and theoretical expectation. This could mislead people to think that the bins are independent, 
which is not true. For this reason, this document does not deal further with unfolded spectra. 



2 Some ways to present deviations in Poisson distributed data 

The goal of this section is to present a logical progression of plotting options, which will lead to the final proposal. We 
examine the different options in detail because they continue to be used despite many people know their drawbacks. 
All the examples will use the observed (from Monte Carlo simulation) and expected event counts of Table 1 and 
assume Poisson-distributed bin populations, which is the most common case in high-energy physics (Sec. 1.3 explains 
the treatment of binomial data) . Two additional features have been introduced in the observed data, an excess and a 
deficit, to make the illustration more interesting. In case the data follow neither a Poisson nor a binomial distribution, 
one needs to compute the p-value according to the probability distribution which the data are assumed to follow. 
Besides the different p-value definition, the rest of this discussion still applies. 

2.1 The D/B ratio 

A simple, and very common practice (e.g.. Fig. 13 in Ref. [6] and Fig. 1 in Ref. [7]), is to plot the ratio of data over 
expectation (D/B), or subtract 1 to plot their relative difference, i.e. {D — B)/B. As shown in Fig. 2, clearly this is 
the worst way to compare two histograms, especially when the bin contents span several orders of magnitude. It is 
possible to hide significant discrepancies, as it will become obvious. 

The biggest disadvantage of this method is that it does not show the statistical significance. It gives the impression 
that low-population bins are fluctuating more significantly than high-population bins, which is not true in this case. 
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observable observable 

Fig. 2. The relative difference between data (D) and expecta- Fig. 3. The quantity (D — B)l\fB approximates the statis- 
tion (_B): [D — B)/B is not a good indicator of the statistical tical significance in bins with large population. However, this 
significance of the observed deviations. approximation fails when there are only a few entries. 

The last non-empty bin is off scale, and if we un-zoom the vertical axis to include it, then the other features will 
become even harder to see. Although the relative difference cannot be smaller than -1 (for D = 0), it has no upper 
limit, therefore, the absolute values of negative and positive deviations are not comparable. Bins with D = would 
all be at relative difference —1, regardless of S, and are typically not shown. 

The only advantage of this method is that excesses of data appear above 0, and deficits below, which is intuitive. 

Although this method is inappropriate to convey the statistical significance of deviations, it might still be used 
when a statistical comparison is not desired. For example, one may wish to show simply that the data lie within 
5% from the expectation, even if a 5% deviation would be statistically significant in regions with many events and 
insignificant in regions with just a few. 

2.2 The {D — B)/^/B approximation 

It is well-known that, for large -B, the Poisson distribution with parameter B is approximated well by a Gaussian 
of mean B and standard deviation ^/~B. This leads to the approximation of the Poisson z- value with , which is 
plotted in Fig. 3. For example, Ref. [8] contains such a plot. 

An advantage of this approach is that we see clearly two significant features: an excess around 200 and a deficit 
around 400. Intuitively, excesses point "up" and deficits "down". Indeed, the quantity shown, in bins with large 
expectation, is a good approximation of the actual Poisson z-value. 

The main disadvantage is that, in low-population bins, is not a good approximation of the true 2;-value. We 

V B 

will see soon that the last bin with _D = 2 is not as significant as it appears in Fig. 3. Similarly, the empty bins in 
the end of the spectrum do not actually have deficits of z-value between and 1, as is indicated by the small negative 
bars. Such z-values would correspond to p-values < 0.5, but, as we are about to show, these bins contain deficits of 
p- value > 0.5, which the approximation in Fig. 3 fails to describe. 

2.3 Plotting signed z-values 

Instead of approximating the z-value with , it would be better to plot the exact z-value in each bin. If a bin 

contains a significant excess of data (with p-value < 0.5), it will have a z-value > which, when plotted, will give 
correctly the impression of an excess. If a bin contains a significant deficit of data (with p- value < 0.5), it will have 
a z-value > as well, so, to avoid giving the impression of an excess in bins with deficits, instead of z-value we plot 
—z-value in these bins. The result is in Fig. 4. 
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observable observable 

Fig. 4. Showing the 2:- value corresponding to the p- value in Fig. 5. The final proposal. Plotting the z-value only for bins 
each bin. In bins with D < B, the — z-value is plotted, to with p- value < 0.5, with positive sign for excesses and negative 
distinguish excesses from deficits. for deficits. 

This is a more accurate representation of significance than , but there is still something confusing in bins with 

low statistics. We know that the last few bins contain deficits (D = 0), however, they appear in Fig. 4 as if they were 
excesses. This happens because they originally had negative z-values, indicating very insignificant deficits, so, with the 
sign-flipping they appear like excesses. The less significant the deficit, the larger the p- value, and the more negative 
the z-value, which gives the impression of an even larger excess when plotting —z-value. The same problem, but in the 
opposite direction, is observed in bin 33, around 800. It contains a very insignificant excess with z-value < 0, and since 
it is an excess this z-value is plotted without changing sign. The result points "down" , giving wrongly the impression 
of a deficit. 



2.4 The final proposal: plot signed z-values only if p-value < 0.5 

In conclusion, our proposal is to plot z-values as in section 2.3, which are negative in bins with data deficits, with the 
exception of bins with p-value > 0.5, whose significance should be simply not shown. The result is shown in Fig. 5. 

The justification for not showing the significance of bins where p- value > 0.5 is that such bins agree perfectly with 
expectation: such deviations are totally uninteresting. In addition, trying to present their z-values could be misleading, 
as explained in section 2.3. 

The advantages of this proposal are the following: 

— The z-value is accurate, unlike the (D — B)/^/~B approximation. 

— Positive values represent excesses of data over the expectation, and negative values correspond to deficits. 

— No feature is hidden if it is worth showing, namely if it has a p-value < 0.5. 

— Bins with high and low statistics (like bin 35 with D = 2 and bin 32 with D — 0) are treated in the same way. 

— If the significance is not shown for a bin, it means that there is no relevant deviation from the expectation. 

— What we propose is easy to implement, using existing ROOT functions as we show in sections 1.1, 1.2 and 1.3. 



3 Including the theoretical uncertainty 

Any theoretical uncertainty in the reference value will affect the significance of the observation. For example, in case of 
a Poisson process with parameter B, a non-zero uncertainty on B will decrease the significance of the difference between 
D and B. Because realistic use cases always involve some (perhaps small) theoretical uncertainty, it is important to 
study its effects on our estimate of the significance of the result. A full treatment would require to solve an inference 
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Fig. 6. The error bars in the inset indicate the range by which 
the significance varies when recomputed after shifting the ex- 
pectation by ±f standard deviation, corresponding to the error 
boxes shown in the main plot. This is not the correct way of 
accounting for theoretical uncertainties. 
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Fig. 7. The inset shows the significance computed by neglect- 
ing the theoretical uncertainty on the expectation (red his- 
togram) and by including it (blue histogram). Including addi- 
tional sources of uncertainty decreases the significance. 



problem. However, here we limit ourselves to the aspects which are relevant when making plots, and consider only the 
total theoretical uncertainty, without addressing all possible components which may contribute to it. 

Of course, in some cases it makes sense not to show the effect of theoretical uncertainties. For example, this can 
be done if the theoretical uncertainty is negligible compared to the statistical uncertainty, or if the deviations without 
theoretical uncertainty are already not significant enough to deserve further investigation. On the other hand, before 
claiming than an excess is really statistical significant one should always check what happens when including the total 
uncertainty on the expectation. 

At first sight, a simple procedure to warn the reader about the effects of the theoretical uncertainty on B would 
seem to plot the significance shown in Fig. 5 with "error bars" . The upper and lower deviations with respect to the 
significance computed above could correspond to the result of recomputing the significance with the parameter B 
shifted by ±1 standard deviation (assuming symmetric uncertainty). An example of this is in Fig. 6, where the relative 
uncertainty in bin i has been set to 10~^ • i^. The problem with this approach is that it misses the fundamental point: 
any additional uncertainty will decrease the significance of the observed deviation. 

This can be shown with a formal approach, which would correspond to a Bayesian treatment of the uncertainties. 
It consists of computing the p-value with the marginal model obtained after integration over the parameter, whose 
prior describes the experimenter's degree of belief on the allowed range of values. This is explained in detail in the 
next two sections for the Poisson and binomial models, and it is illustrated by Fig. 7, which should be compared to 
Fig. 6. 



3.1 Poisson model with uncertain parameter 

In any bin the probability of observing n events is given by the Poisson distribution with parameter y > (the 
expected number of events), which is not certain. One assumes to know its best estimate (which we take to be the 
expectation B of Y) and "uncertainty" (which we take to be the square root S of the variance of Y) . In most cases 
B and S are the result of some auxiliary measurement. We choose to model the uncertainty on the Poisson paramter 
with a function belonging to the Gamma family, which is the conjugate family for the Poisson model. This allows to 
perform the integration analytically. A unique Gamma density corresponds to the pair {B, S). If we write the Gamma 
density as 

Ga(a;|a, b) = x^-^ e"''" (12) 
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the expectation and variance are 

E=y=B and V = ^ = 
b 0"^ 

from which we obtain the Gamma parameters 

a = ^ and b^ — . 

In order to evaluate the probability of observing n events, when the expected yield is described by the Gamma 
density above, we make use of the marginal model which is obtained by integrating over all possible values of Y. The 
result is the Poisson-Gamma mixture (also known as negative binomial) 

b*^ r(n + a) 

P{n\a, b) = Poi(nlr) Ga(r|a, b) dY ^ nlH + b)^+'^ ^^^^ 

(the Gamma function is available in ROOT as r{x) = ROOT: :Math.: :tgaiiima(x)). We now proceed as above, by defining 
the p- value and z-value^ in terms of the marginal model (13): 



p- value = < 



P(n>nobs) = l- Pin\a,b) ^ / 7V(a;; 0, 1) da; , riohs > B 

^ _n ^ 2-valuc 



^'-obs /" 2- value 

P{n<nobs)^YP{n\a,b)^ N{x]Q,l)dx , Uobs < B 



(14) 



(some care needs to be taken when summing over the marginal model (13): it is recommended to make use of the 
recursive relations which allow to compute the n-th term starting from the term n — 1). Finally, we plot the 2;- value 
only if p < 0.5 as explained above, to obtain the results shown in Fig. 7. 

Often the theoretical uncertainty is modeled with a Gaussian distribution, instead of the Gamma density used 
here. Provided that the part of the Gaussian which extends to negative (unphysical) values is negligible, one obtains 
a result (by numerical integration) which is practically identical to what we obtained above (analytically) using 
Gamma. However, in cases where the negative Gaussian tail is non-negligible, the Gaussian needs to be truncated at 
zero. This makes its mean parameter different from its actual mean, and its width parameter different from its actual 
standard deviation. In addition, the truncated Gaussian gives non-null probability to the value of zero (see Ref. [10] for 
additional comments). For these reasons, a truncated Gaussian is a poor choice for the probability model. It is highly 
recommended to use a density function with the correct domain, as the Gamma density, which has the additional 
advantage of making the integration possible analytically. 



3.2 Binomial model with uncertain parameter 

Let the observed ratio be kahs/n in a bin for which we expect a binomial distribution with parameter e £ [0, 1] which is 
uncertain. We choose to model our knowledge about e by means of a Beta distribution with known parameters a, b. The 
Beta functions constitute the conjugate family for the binomial model [9], which allows us to perform the integration 
analytically. The parameters a, b can be found by with the method of moments or numerically. We assume that they 
have been determined in order to best represent the result of some auxiliary measurement. Then the probability of 
observing k events passing the selection out of the initial n events is given by the marginal model 

where the last integral is the definition of the Euler's Beta function B{k + a,n — k + b). Hence 

uni h\ (n\ B{k + a,n-k + b) 

P(fc|n, a, 6) = (^^ j (15) 

In ROOT, the Beta function B{x,y) ~ F (x) F (y) / F {x + y) is available as 

B{x,y) = ROOT: :Math: :beta(x,y) 



^ Defined as in eq. 14 and 16, the z-value is automatically negative for deficits. 
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which makes the marginal model (15) easy to compute, although summing over it is better achieved if the recurrence 
relations which connect P{k\n,a,b) to P{k — l|n,a, 6) are exploited. 

The p-value and z-value are then defined in terms of the marginal model (15) as 



Finally, one plots only the z-values which correspond to p > 0.5 as explained in the previous sections. 

Similarly to the Poisson case, the choice of a Gaussian distribution to model the theoretical uncertainty is problem- 
atic, because in the binomial case the domain is limited to the interval [0, 1], so the Gaussian may have to be truncated 
on both sides. The Beta density is strongly recommended because it has the correct domain and it allows for ana- 
lytic integration. In cases where the Gaussian distribution does not span much below or above 1, it is practically 
indistinguishable from the Beta distribution with the same mean and standard deviation. 



In this article, we propose an improved way of plotting the difference between data and expectation, based on the 
significance of the difference in each bin. We have shown than one can obtain an accurate plot of the statistical 
significance of the deviation of the bin contents from the expectation, obtaining at the same time an intuitive picture 
of the relevant deficits and excesses. This is achieved by computing the exact p- value and, when its value is smaller that 
50% probability, by mapping it into the z-value which gives the deviation in units of Gaussian standard deviations. 
The sign of z-values is always positive for excesses and negative for deficits, making it trivial to decode the result at 
first sight. 

When there are important deviations, whose z-values are above 3 or 4, it is fundamental to check what happens by 
including the total uncertainty on the expectation. This always^ lowers the actual significance and provides a result 
which is closer to what one would have obtained by performing a proper statistical inference. Here the focus is only 
on methods which can improve illustrations. However, it is important to notice that our formulae for the p-values in 
the presence of uncertainty on the expectation coincide with the solutions pr for the Poisson case and pBi for the 
binomial case, computed by Cousins, Linnemann, and Tucker [10] using a completely different approach. These authors 
have performed frequentist evaluations of the performance of such estimators and concluded that they are optimal 
for statistical inference. We repeat that we do not focus on statistical inference here, but this agreement between 
completely independent derivations makes us more confident about the validity of our proposal. 
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